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1 INTRODUCTION 

There is overwhelming evidence that perturbation series in quantum held theory and related disciplines 
diverge (see for example and references therein). Consequently, resummation techniques, which allow 
to associate a hnite value to a divergent series, are needed if divergent perturbation series are to be 
used for numerical purposes. The best known and most often applied resummation techniques are Pade 
approximants Q and the Borel method Recently, another summation method - the so-called delta 

transformation (Eq. (8.4-4) of j^]) - has gained some prominence as a summation method in various 
domains of physics P~|l2"|]. There is evidence |?],|l0| that the delta transformation is able to sum divergent 
series whose coefficients c„ grow essentially like n!, (2n)! and even (3n)!. This is not achievable by 
employing Pade approximants [ jT3| . 

In the case of both Pade approximants and the delta transformation, only the numerical values of a 
finite number of partial sums of the divergent input series are needed as input data. Within the framework 
of the Borel method, it is also necessary to know the large order asymptotics of the coefficients of the 
divergent power series. Consequently, the Borel method is slightly less general than the other two methods. 

Rational approximations to (divergent) power series have an interesting feature: It is possible to 
extrapolate higher-order coefficients that were not used for the construction of the rational approximants. 
This feature, which was apparently hrst observed by Gilewicz Q|, has so far been used quite extensively 
in the case of Pade approximants for divergent perturbation expansions from quantum held theory and 
related disciplines. We refer to the investigations by Elias, Steele and Chishtie et al. Jl5|-pT[|, by Samuel, 
Gardi, Karliner, Ellis et al. ]22| {4]| a nd by Cvetic et al. J4|-|45|]. An analogous approach also works in the 
case of the delta transformation~pl],[l^] . Recursive techniques for the Pade prediction of unknown series 
coefficients were developed recently p™ . 

It is the intention of this article to discuss - augmenting previous investigations [ fT2"[ - how additional 
information on the large-order asymptotics of the perturbative coefficients can be incorporated effectively 
into resummation and prediction schemes. Our emphasis will be on the Borel-Pade method which was 
introduced by Graffi, Grecchi, and Simon p7| . This is a variant of the Borel method that uses Pade 
approximants for performing the analytic continuation of the Borel transformed series to a neighborhood 
of the positive real semiaxis. Further details on the Borel-Pade method as well as on other resummation 
methods can be found in Section |^. 

We would like to mention here the investigations (based on the Borel method) by Fischer |48|-[50f 
and by Caprini and Fischer J^],|5^] regarding the resummation of divergent perturbative expansions in 
quantum held theory. In |5l]] , Caprini and Fischer use asymptotic information (location of the first IR 
and UV renormalon poles) for the construction of conformal mappings in the Borel plane. Here, we 
describe alternative ways in which the asymptotics of the coefficients of the perturbative expansion can 
be utilized to enhance the effectiveness of resummation procedures and perturbative predictions. These 
improvements, which provide an alternative to the methods presented in |5l|| , are not necessarily restricted 
to the first IR and UV renormalons, but can take advantage of the location of all known poles in the 
Borel plane. Regarding asymptotic properties of perturbative coefficients in quantum field theory we also 
refer to fH||]. 

Cvetic and Yu in jl5| consider the resummation of the real part of the (one-loop) QED effective action 
(or vacuum-to vacuum amplitude) in the presence of background electric and magnetic fields, of which 
the exact nonperturbative answer is known [^6). In the presence of an electric field, the QED effective 
action acquires an imaginary part proportional to the pair production amplitude for real electron-positron 
(or lepton-antilepton) pairs. The real part of the effective action does not constitute the full physical 
solution. The full, complex-valued answer requires an integration in the complex plane, for example along 
the special contour introduced in J57| . Only in the case of the pure magnetic field, which we consider in 
this article, the QED effective action is entirely real. For this particular example, the delta transformation 
employed in Jj"2"| or the Borel-Pade Cauchy principal value method used in [ ff5| provide the full, i.e. entirely 
real, not complex physical solution. In view of the above, we again consider here only the resummation 
of the pure magnetic held. The electric held or the electric held combined with a magnetic held should 
be treated along the ideas introduced in |^7j . 

We briefly discuss the relation of resummation and perturbative predictions: There is evidence (see for 
example ^,|57| and references therein) that the resummation of divergent series is ambiguous, especially 
if these series do not fulfill a Carleman condition (see |58| or Theorems XII. 17 and XII. 19 in |5^]). However, 
the prediction of perturbative coefficients apparently does not suffer from such ambiguities. Perturbative 
predictions should be possible |p_ 2[] even in those cases where the perturbation series is evaluated "on the 
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cut" in the complex plane. For the prediction of unknown perturbation series coefficients, the nonanalytic 
contributions, which are responsible for the ambiguities, are irrelevant. Only the analytic part of the 
function, which is represented by a divergent power series, matters. The terms of an otherwise more 
problematic divergent nonalternating series can be predicted just as well as the terms of an alternating 
series. That is to say, the resummation of a divergent series is not necessarily unique (see the different 
integration contours in |p7|), but the perturbative coefficients, which can be extrapolated and predicted via 
rational approximants, are uniquely determined even though the complex integration along the different 
contours in [ p7| leads to different results for the nonperturbative, nonanalytic contributions. 

This paper is organized as follows. In Section |^ we present a short account of various resumma- 
tion methods for divergent series. The exploitation of additional available asymptotic information in 
resummation algorithms is discussed in Section [j[ Asymptotically optimized perturbative predictions 
(i.e., predictions of higher-order unknown perturbative coefficients) are discussed in Sections |]and|^ for 
various example cases. We conclude with a discussion of the results in Section |^. 

2 A REVIEW OF RESUMMATION METHODS 

In this Section, we provide a condensed description of the resummation methods under consideration in 
this article. Let 

oo 

(1) 

be a (formal) power series for some function /. Then, we define the (re, A)-generalized Borel integral 
transform according to 



f(z) = / t A - 1 S^ A )(/;^ K )exp(-i)df (2) 
Jo 

poo 

= z- x ' K s x - 1 B^ x) {f;s K )e^(-s/z 1 / k )As. (3) 



Here, 

oo 

is the (ft, A)-generalized Borel transformed series of the power series (0) for f(z). For re = A = 1, we 
recover the usual formulas for the Borel transformation: 

/•OO 

f(z) = / B(f-zt)exp(-t)dt (5) 
Jo 

1 f°° 

= - / B(f;s)ex V (-s/z)ds, (6) 
z Jo 

OO 

B(f;z) = B^\f-z) = J2^J zU - ( ? ) 

v=0 ' 

There exists an extensive literature on the Borel method in general and on physical applications in 
special. Any attempt to provide something resembling a reasonably complete bibliography would clearly 
be beyond the scope of this article. Let us just mention that recent monographs on the Borel method 
and related topic were published by Shawyer and Watson J60| a nd Sternin and Shatalov Q . 

Let us assume that the coefficients 7„ of the power series (|l|) possess the following large order asymp- 
totics, 

7„ ~ A T(nn + \)B n n -> oo , (8) 

where A, B, n, and A are suitable constants. We say that a (k' , A') -generalized Borel method for some 
power series (Q) is asymptotically optimized if the parameters k! and A' agree with the parameters n and 
A in the large-order asymptotics (j^) for the series coefficient 7„. Thus, the leading (hyper)factorial growth 
of 7„ is exactly canceled out in this case. If we know in addition the parameter B in (||), then we can 
immediately deduce that the asymptotically optimized Borel transformed series (Q) possesses a pole at 
z = l/B. 
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The most difficult computational problem, which normally occurs in the context of a Borel summa- 
tion process, is the construction of an analytic continuation for the Borel transformed scries (|). If the 
coefficients j n of the power series (Q) satisfy (Q), then the (k, A)-generalized Borel transformed series ([IJ) 
has a nonzero but finite radius of convergence. In order to be able to do the integration, we now need 
an analytic continuation which extends B^ K '^'(f;z) from the interior of its circle of convergence to a 
neighborhood which contains the whole positive semiaxis. In this article, we emphasize the Borel-Pade 
method which was introduced by Graffi, Grecchi, and Simon p7[ and which accomplishes the analytic 
continuation by converting the partial sums of the Borel transformed series (Q) to Pade approximants. 
It may happen that the Pade approximants thus constructed exhibit poles along the positive real axis. 
In such a case, the function / is - strictly speaking - not Borel-Pade summable. However, by further 
generalizing the integrals (0) and ([?]) - either via a principal- value prescription (6^,|63|, or by employing 
conformal mappings |3l|,|52f, or via the special integration contours used in f57j[ - it may nevertheless be 
possible to associate a finite value to the divergent power series (Q) for /. 

In the case of Pade approximants we use the notation and conventions of the monograph by Baker 
and Graves-Morris Thus, a Pade approximant [l/m]f(z) to f(z) corresponds to the ratio of two 
polynomials Pi(z) and Q m {z), which are of degrees I and m, respectively, in z: 

r,, i / \ Pi{z) Va+P\z + ...+p t z l 

\l/m\f(z) = — — r— = . (9) 

The polynomials Pi{z) and Q m (z) are constructed so that the Taylor expansion of the Pade approximation 
agrees with the original input series (|l|) up to terms of order I + m in z, 

f(z) - [l/m] f (z) = 0(z i+ ™ +1 ), *-0. (10) 

This asymptotic error estimate leads to a system of linear equations by means of which the coefficients 
Po, Pi, . . . , pi and qi, q2, . . . , q m in (Q) can be computed. However, there are also several algorithms 
which permit a recursive computation of Pade approximants. A discussion of the merits and weaknesses 
of the various computational approaches can for instance be found in Section II. 3 of the book by Cuyt 
and Wuytack @. 

An example of such a recursive algorithm is provided by Wynn's epsilon algorithm ]65| ] : 

- 0, 4 n) =*»» nGN Q , (11) 

e& = 4- + i 1) + V[4 n+1) -4 n) ]. k,n & N Q . (12) 

Wynn |6E| showed that if the input data s n for the epsilon algorithm are the partial sums of the (formal) 
power series ([j]) for some function f(z) according to 

n 

Sn = fn(z) = ^2luZ U , (13) 
v=0 

in) 

then the elements e 2j J with even subscripts are Pade approximants to / according to 



4? - [n + k/k] f (z). (14) 

In contrast, the elements 4fc+i with °dd subscripts are only auxiliary quantities which diverge if the 
whole process converges. 

If one tries to sum a divergent power series or to accelerate the convergence of a slowly convergent 
series by converting its partial sums to Pade approximants, it is usually a good idea to use either diagonal 
Pade approximants, whose numerator and denominator polynomials have equal degrees, or - if this is not 
possible - to use Pade approximants with numerator and denominator polynomials whose degrees differ 
as little as possible. If we use the epsilon algorithm for the computation of the Pade approximants, then 
Eq. (|lj) implies that we then obtain the following staircase sequence in the Pade table (see Eq. (4.3-7) 
of |): 

[0/0], [1/0], [1/1], . . . , [v/v], \v + [v (15) 

This staircase sequence exploits the available information optimally if the partial sums f m (z) with m > 
are computed successively and if after the computation of each new partial sum the element of the epsilon 
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table with the highest possible even transformation order is computed. With the help of the notation [xj 
for the integral part of x, this staircase sequence can be written compactly as follows: 

4£/2? /2l) = [n-[n/2]/In/2]] / (2) ) n = 0,l,2,.... (16) 

The asymptotic error estimate ( |To| ) implies that all series coefficients, which are employed for the 
computation of the Pade approximant [Z/m]/(z), are recovered by a Taylor expansion. Consequently, the 
higher order derivatives of the Pade approximant provide predictions for "unknown" series coefficients, 
i.e. to those series coefficients that were not used for computation of [l/m]f(z). 

There is an enormous amount of literature on Pade approximants in general as well as on their 
application in theoretical physics. Let us just mention that the popularity of Pade approximants in 
theoretical physics can be traced back to a review by Baker |36|, that the monograph by Baker and 
Graves-Morris [Q] is the currently most complete source of information on Pade approximants, and that 
an account of the historical development of Pade approximants and related topics is given in a monograph 
by Brezinski p7| . 

The intense research on Pade approximants during the last decades of course also showed that Pade 
approximants suffer - like all other numerical techniques - from certain limitations and weaknesses. 
For example, Pade approximants arc in principle limited to convergent or divergent power series, but 
cannot help in the case of many other types of slowly convergent or divergent sequences. Moreover, 
Pade approximants are either not useful or cannot be applied at all in the case of power series whose 
coefficients 7„ grow like (2n)! or even (3n)! jt3|. Consequently, the intense research on Pade approximants 
also stimulated research on related techniques, the so-called sequence transformations. 

Let us assume that {s n }n° = o ^ s a sequence, whose elements may for instance be the partial sums of 
an infinite series according to s n — Y) k _Q a^. A sequence transformation is a rule which maps a sequence 
{ s n}%Lo to a new sequence {s' n }^ =0 with hopefully better numerical properties. In this terminology, Pade 
approximants are just a special class of sequence transformations since they transform the partial sums 
of a (formal) power series to a doubly indexed sequence of rational approximants. 

If {s n }n°=o either converge to some limit s as n — > oo or can be summed to the generalized limit s in 
the case of divergence, then a sequence element s n can for all n > be partitioned into the (generalized) 
limit s and a remainder r n according to 

s n = s + r n . (17) 

Normally, a sequence transformation will not be able to determine the (generalized) limit s of {s n }J° =0 
exactly. Thus, the elements of the transformed sequence {s^}^o can a l so be partitioned into the (gener- 
alized) limit s and a transformed remainder r' n according to 

S'n = « + < , (18) 

and the transformed remainders will in general be different from zero for all finite values of n. 

In the literature on convergence acceleration it is said that a sequence transformation accelerates 
convergence if the transformed remainders {r^}^_ vanish more rapidly than the original remainders 
{r„}^L according to 

lim T JSl = ii m £n_l£ = o, (19) 

rwoo r n ra— >oo s n — S 

and a divergent sequences, whose remainders r n do not vanish asn-*oo, is summed to its generalized 
limit s if the transformed remainders r' n vanish as n — > oo. 

Thus, a sequence transformation essentially tries to eliminate the remainders r n from the sequence 
elements s n as effectively as possible. Since, however, an in principle unlimited variety of different remain- 
ders can occur, it is necessary to make some assumptions - either explicitly or implicitly - which provide 
the basis for the construction of a sequence transformation. A detailed discussion of the construction of 
sequence transformations as well as many examples can be found in the book by Brezinski and Redivo 
Zaglia |H or in §. 

Normally, the assumptions being made are incorporated into the transformation scheme via model 
sequences, whose remainders possess a particular simple structure and can be expressed by a finite number 
of terms: 

k-i 

s n = s + } j Cj<pj(n). (20) 

k=0 
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Here, the Cj are unspecified coefficients, and the ipj (n) are by assumption known functions of n. 

The elements of this model sequence contain k + 1 unknown, the (generalized) limit s and the k 
coefficients Cj with < j < k— 1. Since all unknowns occur linearly, it is possible to construct a sequence 
transformation T - if necessary via Cramer's rule - which is exact for the elements of this model sequence 
according to 

T = T (s„, s„ + i, . . . , s n +k, ) = s, (21) 

if applied to the numerical values of k + 1 sequence elements s n , s n +i, ■ • • j •?„+?;. 

Of course, simple model sequences of that kind normally do not occur in practical problems. However, 
their elements provide at least for sufficiently large values of k reasonably accurate approximations to the 
elements of the more realistic sequence 

oo 

S n = S + ^ Cjfji 11 )- ( 22 ) 
3=0 

If we now apply this sequence transformation T to the numerical values of k + 1 elements of the sequence 
(p2|), then we have no reason to assume that T might produce its exact (generalized) limit s. However, 
a more detailed mathematical analysis of the transformation process normally reveals that T eliminates 
the first k terms Cjipj(n) with < j < k — 1. Thus, the transformed remainder r' n starts with tpk(ri) 
instead of <po(n), which for sufficiently large values of k normally constitutes a significant improvement. 

Most sequence transformations can be constructed on the basis of model sequences of the type of Eq. 
(f20|). For example, Wynn fj5l could show that his epsilon algorithm is exact for model sequences whose 
remainders can be expressed as a linear combination of exponential terms according to 

fe-i 

S n = s + ~ c i A " • ( 23 ) 

3=0 

Concerning the Xj it is only assumed that they are different from zero and one and ordered according to 
magnitude, i.e., Xj ^ 0, 1 and |Ao| > |Ai| > |A&_i| > 0. Thus, if the numerical values of 2k + 1 elements 
s n , s n +i, ■ ■ • , s n +2k of this model sequence are available, then the epsilon algorithm is exact according 
to 

4? = S. (24) 



Moreover, Wynn constructed in Theorems 16 and 17 of |69| asymptotic expansions (n — * oo) for the 

>) 

•2k 



transformed remainders r' n created by the application of to the elements of the sequence 



z 3 A" , (25) 



S n = S + 2J 

3=0 

which is an obvious generalization of the model sequence ( |23| ) . He showed that the transformed remainders 
r' n are proportional to A" which corresponds to an elimination of the first k exponential terms CjA™ on 
the right-hand side of (f25|). Since the Aj are by assumption ordered in magnitude, this constitutes a 
significant achievement. Consequently, Wynn's epsilon algorithm is asymptotically optimal for sequences 
of the type of (^5|). This means that no other sequence transformation, which also uses only the numerical 
values of the elements of the sequence ( p5| ) as input data, can produce a better asymptotic (n — > oo) 
truncation error. 

Levin ]70[] introduced a class of sequence transformations which are exact for model sequences of the 
following type: 

S n = S + U) n Z n - (26) 



Here, u> n is an estimate for the truncation error f n , and z n is a correction term. Levin |70| assumed that 
z n can be expressed as a truncated power series in l/(n + where £ is a positive shift parameter. In 
Sections 7 - 9 of ||, several other sequence transformations were constructed which are also exact for 
the model sequence ( p6[ ) but make different assumptions about the correction terms z n . The remainder 
estimates lo„ introduce additional degrees of freedom in the construction of the sequence transformation as 
compared to Pade approximants. One may draw an analogy between sequence transformations and Pade 
approximants on the one hand and the Gaussian integration and the Simpson rule on the other hand; 
the variable integration nodes and weight factors of the Gaussian integration yield additional degrees 
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of freedom which may be used in order to construct a potentially much more powerful algorithm for 
numerical integration. 

In the following text, we will concentrate on sequence transformations which assume that z n can be 
expressed as a truncated factorial series (Section 8 of [|): 

fc-i 

3=0 

Here, (n + ()j = T(n + C + j)/P( n + C) is a Pochhammer symbol, and £ is a positive shift parameter. 
The assumption (ETj) implies that the sequence transformations derived in this way are particularly well 
suited for sequences satisfying 

oo 

s„ = s + u> n Cj/ (C + ri)j . (28) 

3=0 

It is a priori not obvious that the ratio [s n — s]/u> n = r n /u> n can be expressed as a factorial series. Nev- 
ertheless, this assumption leads to powerful sequence transformations which are apparently particularly 
well suited for the summation of factorially divergent series [|| [l^]. Let s n = J2k=o °" k ^ e * ne P ar tial 
sums of an infinite series. If the strictly alternate in sign and decrease monotonously in magnitude, 
then the best simple estimate for the truncation error r„ = — Y^k= n +i a k 1S the nrs t term a n +i neglected 
in the partial sum s n . Moreover, the first term neglected is also the best simple remainder estimate for 
many factorially divergent alternating series (see for example Theorem 13-2 of [61). The mathematical 
structure of a factorially divergent series is expected for the perturbative expansions in quantum field 
theory. Further arguments supporting the general applicability of the delta transformation to the series 
of the mathematical structure as expected for quantum field theory will be discussed in |86fl . 

If we now combine the assumption that z n should be a truncated factorial series according to ( f27j ) 
with the remainder estimate 

uj n = As n = a n+ i , (29) 

which was introduced by Smith and Ford fnfl , we obtain the delta transformation which is defined by 
the following ratio of finite sums (Eq. (8.4-4) of ||): 



A fc [(C + n)fc_ia n /As n ] 
A k [(C + n) fe _i/As„] 



3=0 



k\ (C + n + j)fe-i s n+J 



3=0 



jj (C + n + k)k-i As n+j 

. (30) 



jj (C + n + fc)fe_i As n+j 



Here, the same notation as in (|] is used. Thus, A stands for the difference operator defined by Ag(n) — 
g(n + 1) — g(n), (a) n = T(a + n)/T(a) is a Pochhammer symbol, fc and n are nonnegative integers, and £ 



is a shift parameter which has to be positive to allow n — in Eq. (30). The most obvious choice, which 
is always used in this article, is £ = 1. 

In Section 8.3 of a simple recursive scheme is described which permits - depending upon the initial 
values - the recursive calculation of either the numerator or the denominator sum of 5f, (£, s n ) . 

In the context of quantum field theory and related disciplines, the delta transformation (pG) may be 
used for the summation of divergent perturbation expansions which are power series in some coupling 
constant. Thus, if we replace the input data s n in (^30|) by the partial sums f n (z) — 'YT v =<ilvZ 1 ' of the 
(formal) power series (Q) for f(z), we obtain a rational expression, whose numerator and denominator 
polynomials are of degrees k + n and k, respectively, in z: 

J2 (-1)' ^ (C + "+j)fc-l **~ J '/n+j (z) 



k 

E (-^ 

3=0 



k\ [C + n+j)k-i z k -i 



4 n) (C=/nW) = ^~~~~k — T : 7—- (31) 



j J (C + n + k) k -! 7„+ 



3+1 
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If the coefficients j n of the power series for f(z) are all different from zero, the rational function ( plf ) 
satisfies the asymptotic error estimate (Eq. (4.29) of JTo|] ) 

f(z) - 4' l) (C,/„(z)) = 0(z k+n + 2 ), z^O. (32) 

This estimate, which is formally very similar to the analogous estimate ( ^0| ) for Pade approximants, implies 
that all terms of the formal power series, which are used for construction of the rational approximant 

(Ci fn(z)) , are reproduced exactly by a Taylor expansion around z — 0. Moreover, the higher order 
derivatives provide just like in the Pade case predictions for those coefficients y n +k+2> 7n+fc+3, ■ ■ • > that 
were not used for the construction of the rational function. 

As already discussed, the power of the delta transformation or of other Levin-type transformations 
results from the fact that an explicit estimate for the truncation error is incorporated into the transforma- 
tion scheme. The truncation error estimate used by the delta transformation is the first term ^ n +iz n+1 
neglected in the partial sum f„(z) — X)"=o T^ 1 " '■ Consequently, for a proper application of the delta 
transformation all coefficients 7„ of the power series for / with n > 1 have to be different from zero be- 
cause otherwise the estimate for the truncation error makes no sense. This restriction also follows directly 
from the ratio (|3f|), where undefined expressions occur if coefficients 7„ with n> 1 are zero (cf. Eq. (II) 
in |45|] which entails divisions by zero). 



3 ASYMPTOTICALLY OPTIMIZED RESUMMATION 

The problem of the resummation of divergent perturbative expansions in quantum field theory and related 
disciplines has been discussed in a number of recent publications, for example in [||[|5^,|5^,^,^3|,^4). We 
investigate here asymptotically optimized resummation methods, i.e. methods which utilize information 
about large-order asymptotics of perturbative coefficients with the intention of enhancing the rate of 
convergence of the resummation algorithm. 

We discuss here possible improvements of the Borel-Pade method on the basis of potentially available 
information about the large-order asymptotics of perturbative coefficients. We concentrate on the partic- 
ular model example discussed recently by Dunne and Hall j55| , by Cvetic and Yu in p5| and by ourselves 
in p^] . We discuss the QED effective action Sb in the presence of a constant background magnetic field. 
The exact nonperturbative result for Sb can be expressed as a proper-time integral: 

oo 

e 2 B 2 f ds f , 1 si / m 2 c \ 
o 

Here, B is the magnetic field strength, and m c and e are the mass and the charge of the electron, 
respectively (this result is given for example in Eq. (4-123) in J75|]). 

The integral representation ( |33| ) for 5b can be expressed as a strictly alternating perturbation series 
in the effective coupling coupling — e 2 B 2 /m^: 

2 e 2 B 2 °° 

Sb = o - 9b 22 Cn 9b ) ( 34 ) 

n=0 

where 

(-l)»+i 4" |g 2n+4 | 
n (2n + 4)(2n + 3)(2n + 2) ' V ' 

Here, 6271+4 is a Bernoulli number. Thus, the perturbation expansion ( |34| ) has the remarkable feature 
that an unlimited number of series coefficients c n arc known analytically. Consequently, this series is 
particularly well suited as a model system for studying resummation methods. 

Next, we utilize the fact that a Bernoulli number with even index can be expressed by a Riemann 
zeta function according to (see Section 23.2. on p. 807 of J76[) 
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Inserting this into and (B5J) yields 



*b = -S^E(-ir(2»+l)!^S^^ (37) 



71=0 

e 2 B 2 



5B S B , (38) 

OVl " 

where in the last line we define implicitly the scaled function S B which is also considered, e.g., in Table 
2 of [i2| . It is a direct consequence of the Dirichlet series (Section 23.2. of |[7q| ) 

oo 

C(s) = ^(m + l)-% (39) 

m=0 

that we have the inequality 1 < ((2n + 4) < ((4) for all nonnegative integers n. Consequently, the zeta 
function does not contribute to the factorial divergence of the perturbation series (34). Thus, the factorial 
(2n +1)! on the right-hand side of ( |37| ) implies that the perturbation series for S B diverges for every 
coupling <7b 7^ 0. Furthermore, it is clear from the representation (B7[) that an asymptotically optimized 
(k, A)-generalized Borel resummation scheme for Sb according to (2^-(^) requires the parameter setting 
K = A = 2. 

We now discuss the construction of the asymptotically optimized Borel transform explicitly. We start 
from the scaled series, 



S b (9b) = £( 16c »)^ = £(-!)" (2n+l)! C ^ " + > g B . (40) 



2C(2n + 4) ^ 

n=0 n=0 

The (2, 2)-generalized Borel transformed series of S B is given by 



^(s B ;z) = £ T ^ 5 § = E(-ir^^^- (4i) 

n=0 \ n=0 

This series can be brought into a form which clearly shows the singularity structure of the Borel trans- 
formed series. For that purpose, we replace the Riemann zeta function by its Dirichlet series according 
to (p9|) and interchange the order of the two infinite nested summations: 



,-.00 oo 

B^(S B ;z) = -1 £ (m + 1) -4£ j_ z/[7r(m + 1)] 2j\ (42) 



2 

m=0 

Thus, 2?( 2 ' 2 ) (S' B ; z) is essentially a superposition of geometric series with arguments z/ [7r(m + l)] 2 . If we 
now use J2kL {—x) k = 1/(1 + x), we obtain: 

B^{S' B -z) = A V (m+1) ' 4 2 . (43) 

^ 4 i + ^/K^ + i)] 2 

This representation shows that the poles of the Borel transformed series £i( 2 ' 2 ) (S* B ; z) are located along 
the negative real axis according to 

z = - n 2 ir 2 , (44) 

where n is a nonzero positive integer. Moreover, we obtain the following representations for the QED 
effective action Sb as a (2, 2) -generalized Borel integral accor ding to Eqs. (§) and (§): 

2 R 2 roc 

Sb = -inrSB / i8M(^; 5B t 2 ) exp(-i)di (45) 

2 d2 roc 

= -ITT- / 5 ^ ( 5 b5 ^) exp(- S / 5B /2 ) d S . (46) 
57r Jo 

Cvetic and Yu in jl5) use a (2, 2) -generalized Borel transformed series for the QED effective action, 
constructed according to Eq. (45). As explained in Section |l|, this transformation is asymptotically opti- 
mized in the sense that the leading factorial growth of the perturbative coefficients in Eq. ( |37| ) is divided 
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out. It could appear from the Eqs. (4) and (8) in that a (1, l)-generalized Borel transform is used 
where the nth perturbative coefficient is divided by a factor of nl — T(n + 1). This is, however, not the 
case. Note that in p5[ , the perturbative expansions are written in a very peculiar parameter, which is 

b = eB/m 2 c . (47) 

In normal QED terminology, this would correspond to an expansion in powers of ^/a. By consequence, 
all even-order perturbative coefficients vanish in the analysis presented in Ref. p5| . In this context, one 
may to note that the expression for the delta transformation according to Eq. (11) in |^] is actually 
undefined since it involves divisions by zero. The expansion parameter used here is gB = e 2 B 2 jm\ = b 2 ; 
this parameter is also used in |l^] and in |5q| . 

Due to the special, mathematically compact form of the perturbative coefficients in Eq. ([35j), the 
asymptotically optimized Borel transform (|42| ) is simply a superposition of geometric series. Such a simple 
mathematical structure cannot be expected to be of general importance concerning series occurring in 
quantum field theory. Note that this series can be brought in a form which clearly shows that Wynn's 



epsilon algorithm, which computes Pade approximants according to (14), is optimal, as discussed in 
Section ||. For that purpose, we rewrite the ?ith partial sum of the series (|43|) as follows: 

i^—O u—n+l 

On substituting the Dirichlet series ( |§9| ) into the infinite series on the right-hand side and interchanging 
the order of summations, we obtain 

Thus, the partial sum of the Borel transformed series ( fl3|) possesses the following general structure: 

s n = a + (-1)" c > A " ■ ( 5 °) 

3=0 

This sequence is obviously a special case of the sequence (p5|), for which Wynn's epsilon algorithm is - 
as discussed in Section || - asymptotically optimal. 

The conclusions drawn by Cvetic and Yu in EEI appear to be restricted, at least in part, to the partic- 
ular model example studied in their paper. In this context it should be emphasized that the superiority of 
the delta transformation over Pade approximants if applied directly to factorially divergent series, cannot 
be assumed to persist after the Borel transformation by which the leading factorial divergence is divided 
out. I.e., the delta transformation is more powerful than the Pade technique for factorially divergent 
series, but this finding by no means allows us to conclude that, or in fact has any connection to the 
assumption that the combined Borel-delta technique should be numerically superior to the Borel-Pade 
method. This consideration is relevant for the interpretation of the conclusions drawn by Cvetic and Yu 
with regard to the variant of the Borel method proposed in []45f , which uses the delta transformation for 
the analytic continuation of the Borel transformed series and which is called the Borel- Weniger method 
by the authors of Kf|. 

We return now to the discussion of further improvements of the asymptotically optimized Borcl- 
Pade method. The leading asymptotics do not only permit to modify (optimize) the Borel transform 
accordingly, but indeed it is the leading large-order asymptotics which determine the location of the 
poles in the Borel plane. In view of Eq. (pf3|), the singularities of the function £>( 2 ' 2 ) (S B ;z) defined in 
( p4| ) are at z = —n 2 ir 2 , i.e. along the negative real axis. This is where one would expect them to lie in 
(distant) analogy to the renormalon theory ff% . 

Using the information on the location of the poles, it is possible to construct further improved Pade 
approximants. To this end, we utilize the known location of the poles in order to construct improved Pade 
approximants to the function 23( 2 > 2 ) (S B ;z). Normally, the Borel integral ( [43| ) would be evaluated with 
the upper- or lower-diagonal Pade approximants to B^ 2 ' 2 ^ (S B ; z) in the integrand. We use upper-diagonal 
Pade approximants here, as they can be computed by Wynn's epsilon algorithm according to (|lj). We 
denote by V n (z) the upper-diagonal Pade approximant, 

V n {z) = [[(n+l)/2]/[n/2]] B(a , a)( (z). (51) 
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Table 1: Resummation of the divergent series given in Eq. (^) with g B = 10. 
The numerical data is presented normalized to a number in the interval (0, 1) 
via multiplication by a factor of 100. 



rrl ql q-lll ql 

1 °B,n 1 °B,n 

asymptotically optimized, improved transforms 
see Eq. (||) def. in Eq. © 






-2 


.222 


222 


222 


222 


222 


222 


-2 


.222 


222 


222 


222 


222 


222 


1 


-0 


.779 


860 


343 


938 


511 


943 


-0 


.846 


447 


993 


882 


134 


544 


2 


-0 


.832 


545 


950 


383 


972 


556 


-0 


.807 


698 


096 


764 


310 


129 


3 


-0 


.804 


166 


791 


460 


607 


115 


-0 


.805 


669 


649 


913 


560 


215 


4 


-0 


.806 


776 


251 


699 


410 


322 


-0 


.805 


634 


754 


493 


393 


579 


■5 


-0 


.805 


531 


742 


010 


943 


471 


-0 


.805 


634 


061 


009 


148 


890 


G 


-0 


.805 


700 


473 


754 


628 


870 


-0 


.805 


633 


961 


318 


348 


380 


7 


-0 


.805 


626 


062 


286 


745 


674 


-0 


.805 


633 


975 


558 


025 


628 


8 


-0 


.805 


638 


540 


560 


183 


781 


-0 


.805 


633 


975 


330 


131 


982 


9 


-0 


.805 


633 


322 


312 


161 


338 


-0 


.805 


633 


975 


322 


121 


382 


10 


-0 


.805 


634 


321 


402 


303 


670 


-0 


.805 


633 


975 


321 


669 


649 


11 


-0 


.805 


633 


919 


056 


749 


287 


-0 


.805 


633 


975 


321 


697 


521 


12 


-0 


.805 


634 


003 


326 


479 


025 


-0 


.805 


633 


975 


321 


696 


356 


13 


-0 


.805 


633 


970 


320 


418 


441 


-0 


.805 


633 


975 


321 


696 


294 


14 


-0 


.805 


633 


977 


694 


044 


469 


-0 


.805 


633 


975 


321 


696 


292 



exact -0.805 633 975 321 696 292 -0.805 633 975 321 696 292 



In the upper-diagonal case, we evaluate the transforms TS' B where 

/■oo 

TS B „ = / t V n (g B t 2 ) exp(-i) dt , (52) 
Jo 

and observe numerical convergence of the transform at large transformation order n. When the location 
of the poles is known, we may improve the convergence of the transforms by the following replacement, 

T>n{z)^V' n{z) = ln/n Q ^ Z) , (53) 

J] (l + z/(n 2 ^ 2 )) 

i=l 

where Q n (z) is the [[[(n + l)/2j/0]-Pade approximant to the function lZ n (z), 

Q n (z) = [i(n + l)/2yO] n Jz), (54) 

and lZ n (z) is given by 

[n/2] 

ftn(*) = II ( X + Z /(" V )) Bi2,2) ( 5 b; z) . (55) 

i=l 

The asymptotic enhancement is possible only if additional asymptotic information is available on the 
perturbative coefficients. Such information may be available (renormalon poles) , but this is not necessarily 
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provided. In this context it should be noted that there is currently no general proof of the assumption 
that the renormalon poles are the only relevant poles in the Borel plan e f79||, but the factorial divergence 
of the perturbative coefficients is a commonly accepted assumption p],|77|,|80|"p5[. 

The pole-structure improved transforms are obtained from (p2|) by the replacement (|5J 



/>oo 

T'Sb.u = / t V' n (g B t 2 ) exp(-t)dt. 
Jo 



(56) 



Similar improvement of the convergence of transforms can also be expected in those cases where the final 
evaluation of the Borel integral proceeds in the complex plane along the integration contours introduced 
m §?]]. 

Further improvement of the rate of convergence is possible by taking the transforms T' S'-q n as input 
data to the epsilon algorithm (|ll]) in order to accelerate the convergence of the sequence of the pole- 
structure improved transforms {T'S^ nl^Lg. The application of the epsilon algorithm defined in Eq. ( pd[ ) 
to the pole-structure improved transforms results in a sequence of upper-diagonal Pade approximants 
which we denote by 

T"Q' _ ,(™-2[[«/2l) ( r 7 \ 
1 D B : n — e 2[n/2]] ' \° ' > 

As input data for the epsilon algorithm, we use 

s n = T'S' B n . (58) 

In a second epsilon transformation we may in turn employ the T" S' B as input data for a further 
application of the epsilon algorithm, 

T m q' — ,(™- 2 [[™/ 2 l) fcQ\ 

1 °B,n — e 2[n/2] ' V oy J 

where we use s n — T" S' B n as input data. This results in a sequence of pole-structure and doubly epsilon- 
improvcd transforms T" l S' Qn . The application of the epsilon algorithm further enhances the rate of 
convergence of the pole-structure improved transforms. 

In Table [l] we present numerical data for the Borel transforms calculated according to Eq. ([32]) (in 
passing we note that these correspond the method proposed in p5[) and the transforms (|59|). The first 
15 transforms calculated according to Eq. ( |52| ) exhibit convergence to 8 significant digits, whereas the 
pole-structure and epsilon improved transforms coincide with the exact result to within 18 significant 
digits. 

The delta transformation is a general-purpose transformation which has been proven to be applicable 
to a wide variety of alternating factorially divergent series [^|-|[^0 12 1. In the context of the delta trans- 



formation, additional asymptotic information could be used in order to modify the remainder estimates 
uj„ defined in Eq. @ (see also Eqs. (7.3-8), (8.2-7), (8.4-1) and (8.4-4) of §). Also, we note a rescaling 
of the perturbative coefficients as a potential source for further improvements . Work along these lines 
is currently in progress and will be presented elsewhere [p6|. 



4 ASYMPTOTICALLY OPTIMIZED PREDICTIONS 

We refer here to the predictions of unknown higher-order perturbative coefficients as perturbative predic- 
tions or perturbative extrapolations. As outlined in Section [| and |2q] , these extrapolations are obtained 
by reexpanding certain rational approximants in powers of the coupling parameter. The next higher- 
order term obtained after the reexpansion can then be interpreted as a prediction for that perturbative 
coefficient. The rational approximants discussed in Section ^| fulfill accuracy-through-order relations, i.e., 
upon reexpansion in the coupling parameter, all the perturbative terms used for the construction of the 

) for Pade approximants and Eq. ( p2] ) for the delta 
transformation] . The coefficients of the Borel transformed series are related to those of the input series 
by Eq. (g). Therefore, we can either predict the perturbative coefficients of the original series or the 
coefficients of the Borel transformed series. 

We consider here the asymptotic improvement of three different prediction methods: (i) asymptotically 
optimized predictions based on the combination of Borel and Pade techniques, (ii) predictions based on 
the delta transformation, and (iii) the direct application of Pade approximants to the perturbation series. 
We consider here the following improvements beyond reexpansion of the rational approximant, 



rational approximant are reproduced [see Eq. (10 
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1. A-postcriori corrections. These are further corrections to the perturbative predictions obtained by 
estimating not only the coefficient, but also the probable error in making that estimate. Similar 
methods in the context of Pade approximants have been introduced, e.g., in p!q ]. 

2. Fixing poles. In the context of the asymptotically improved Borel-Pade method, it is possible 
to improve predictions if the leading and subleading large-order asymptotics of the perturbative 
predictions are known. These asymptotics determine the location of the poles in the Borel plane, 
which can be put in by hand (see also Eqs. (p3|)-([55|) in Section ^ and, in part, |92j|). 

3. Renormalization group. The renormalization group can be used to enhance perturbative predictions 
for certain classes of diagrams; this has been used e.g. in Q for the anomalous magnetic moment 
of the muon. 

It is natural to assume that combinations of these techniques should be investigated where appropriate. 

The basic idea of a-posteriori corrections is as follows. The errors made in the "prediction" of lower- 
order coefficients are available by the time we come to higher order, so they may be utilized for an estimate 
of the error which is to be expected in a prediction of the next higher-order coefficient. We denote the 
predictions which are obtained by extrapolating the coefficients and the "prediction errors" (in contrast 
to the coefficients alone) by the term a-posteriori improved predictions because the further correction due 
to the extrapolated error is applied after the reexpansion of the rational approximant which yields the 
"first-order" prediction. The a-posteriori improvement of predictions is useful in both lower and higher 
orders of perturbation theory. In higher orders, the transient, pre-asymptotic behavior of the perturbative 
coefficients has died away, and the extrapolations of the coefficients as well as the a-posterori corrections 
become more accurate. A particular merit of the a-posteriori corrections is the fact that they can be 
applied to any of the prediction algorithms proposed above, in order to achieve an improvement of the 
prediction beyond the reexpansion of the rational approximant used. 

We will be investigating in the sequel the Borel transform of the QED effective action defined in 
( p3| ) . In order to investigate the extrapolation of coefficients of the Borel transform, we define auxiliary 
quantities (coefficients) C2 n +i by the relation 

16 On (cn , 
(2n + 1)! 

where the c„ are defined in Eq. (|35|). We additionally set c<x n = for even-order coefficients. In terms of 
the coefficients Cj(p) which are defined in Eq. (5) in ifjsf , the C2 n +i are given by C2 n +i = (— l) n C2n+i(0). 
The coefficient C2 n +i, written in terms of the Bernoulli numbers, reads [see also Eq. (p5|)] 

» _ , 2 2 "+ 4 |ff 2 » +4 | , , 

C2 " +1 - ( ~ 1} (2n + 4)! ■ (61) 



We concentrate here on the coefficient C13 defined in Eq. (|6G|) and we discuss how the prediction 
for the coefficient £13 can be improved on the basis of a-posteriori corrections and other asymptotic 
improvements. We define correction factors £ n by 

On £n On (62) 

where c„ is the exact nth order coefficient and c„ is the estimate obtained by reexpanding the rational 
approximant which is used for the prediction. In the case of Borel-Pade approximants, these would be 
Pade approximants applied to the Borel transform of the QED effective action (|4^). Specifically, for the 
prediction of the nth perturbative coefficient, these are the approximants [(n— l)/2]/[n/2] for the lower- 
diagonal and [n/2]/[(n— l)/2] for the upper-diagonal case. For the prediction of C13, the previous errors 
made in the "prediction" of 67, eg and en can be analyzed. Note that the exact values of £7, eg and en 
must be assumed as available, exploitable information by the time we try to predict £13. 

An estimate for £13 can be obtained for example by fitting the natural logarithms of the quantities 
£7 — 1, £9 — 1 and £n — 1 with a linear model in order to obtain an estimate for £13 — 1. This linear fit 
of the logarithms of the £j is based on the empirical observation that relative errors of the predictions 
decrease exponentially in higher order, a phenomenon which has been observed in a number of applica- 
tions, including variants of anharmonic oscillators. In the context of Pade approximants, a similar error 
dependence has been conjectured (see The leading coefficient of the decay of the relative errors may 
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Table 2: Prediction of perturbative coefficients £13 defined in Eq. (|60|). 



method result relative error 

Asympt. opt. [5/6] -Borel-Pade 2.221 459 447 36 x 1CT 8 6 x 1(T 7 
(see P5| ) 

Asympt. opt. [6/5]-Borel-Pade 2.221 454 724 11 x 1CT 8 3 x 10~ 6 
(for comparison) 

Asympt. opt. [6/5]-Borel-Pade 2.221 460 221 71 x 1CT 8 3 x 1(T 7 
with a-posteriori correction 
(this work) 

Asympt. opt. [6/5]-Borel-Pade 2.221 460 950 35 x 10~ 8 3 x 10~ 8 
with one pole (this work) 

Asympt. opt. [6/5]-Borel-Pade 2.221 460 901 25 x 10~ 8 1 x 10~ 8 

with three poles (this work) 

Asympt. opt. [6/5]-Borcl-Pade 2.221 460 880 27 x 10~ 8 5 x lO" 10 

with five poles (this work) 

Asympt. opt. [6/5]-Borel-Pade 2.221 460 878 35 x 10~ 8 3 x 10~ 10 

with five poles + a-posteriori (this work) 

ci 3 : exact value 2.221 460 878 99 x 10~ 8 



depend on the problem considered and on the extrapolation scheme used, but the exponential improve- 
ment of perturbative predictions in higher order appears to be a rather general feature. Details on this 
point will be presented elsewhere j8f|. 

Using a linear least-squares fit of the respective logarithms ln(^7 — 1), ln(^9 — 1) and ln(£n — 1), we 
obtain an estimate of £ 13 — 1 = 2.47 x 10 -6 in the case of upper-diagonal Borel-Pade approximants. This 
leads to the data presented in Table ||. Note that even the crude linear model for the ln(£, — 1) used here 
already doubles the accuracy of the prediction of the coefficient C13 as compared to the plain Borel-Pade 
prediction used, for example, in Note also that an averaging of upper-and lower-diagonal Borel-Pade 
approximants does not improve the situation in favor of the plain Borel-Pade extrapolation. Further 
improvements of the a-posteriori corrected predictions is possible with more elaborate extrapolation 
schemes @. 

We would like to mention that knowledge of the leading asymptotics of the perturbative coefficients is 
required for the construction of an asymptotically optimized Borel-Pade transformation; this information 
is not available for many of the phenomenologically interesting series currently investigated jl5]-^l] . It 
is helpful to note that the delta transformation is (like Pade) a rather general-purpose method for the 
prediction of perturbative coefficients, and that knowledge of large-order asymptotics is not required 
for its construction or application in a particular case. This property is helpful especially in cases of 
practical interest where little is rigorously known about the large-order asymptotics of the perturbative 
coefficients, and where only a limited number of perturbative coefficients are available. A number of 
practically interesting examples of delta-based predictions were discussed in ]l2| , and it was observed 
that the delta transformation yields more accurate predictions than the Pade technique in many cases. 

We now turn to a discussion of topologically new effects in higher orders of perturbation theory and 
improvements of predictions based on the renormalization group. Topologically new effects have caused 
problems for perturbative predictions in the past. We refer to the quartic Casimirs in the QCD beta 
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function |3(], [h^, |8{| and to light-by-light scattering graphs in the tenth order anomalous magnetic mo- 
ment of the muon |mJ. Analogous considerations might hold for the perturbation series investigated 
in p2| ]. The topologically new effects cannot be taken into account by straight extrapolations, nor by 
renormalization-group improved extrapolations, which lead to resummation of certain classes of dia- 
grams. For the anomalous magnetic moment of the muon discussed in pffl, the contribution of the 
topologically new light-by-light scattering diagrams originally analyzed in pl[ could not be reproduced 
by renormalization-group techniques introduced in [ |90| . By reexpansion of the delta approximant to the 
perturbation series for the muon anomaly, an estimate of a^ 10 ^ =711 has been obtained ]87| . If an a- 
posteriori correction based on a combination of the delta transformation and Pade approximants is added 
to this prediction, then the estimate for the 10th order coefficient changes to aj, 10 ' ~ 970 J86|. This im- 
proved estimate is in excellent agreement with the analytically obtained approximate result of 930(170) 
from ]9l| ] , comprising the topogically new effects which are present in five-loop order. 

Now it is of course not permissible to conclude that topologically new effects can be taken into account 
in the general case by considering a-posteriori corrections, and/or that a-posteriori corrections necessarily 
give an accurate estimate for the size of these problematic, topologically new effects. On the other hand, 
the reverse proposition, which is that perturbative predictions are scientifically unsound because of their 
general inability to include topologically new effects, does not appear to be generally valid, either. 



5 THE ANOMALOUS DIMENSION 

We consider the divergent perturbation series for the 7 function (anomalous dimension) of the Yukawa 
coupling as studied in J78|. Specifically, we consider the resummation of the perturbation series for the 
anomalous dimension of a fermion field with a Yukawa interaction gip a ip at d c = 4, which is given in 
Eq. (17) in ]78| ]. This calculation comprises an evaluation of the contribution of all nested self-energy 
diagrams to the anomalous dimension 7 function of the Yukawa theory up to the 30-loop level (an 
analogous analysis is performed in |78| for the (4— e)-dimensional (f> 4 theory). We restrict the discussion 
here to the Yukawa case. With the convention 

the result for the anomalous dimension 7 function as considered in |7§| ] reads, 

00 x, 

WW~EH)"^T fl "- ( 64 ) 

71=1 

The perturbative coefficients G n are listed in Table |[ The coefficients grow factorially in absolute mag- 
nitude; in |78| little change is observed in the quantities 

Sn = 2"-iI> + l/2) (65) 

for large n. The evaluation |7§f| confirms in a concrete, 30-loop calculation the assumption originally put 
forth by Dyson |8S|| that the convergence radius of the quantum field theoretic perturbative expansion is 
zero. For a large number of quantum field theoretic observablcs like the anomalous magnetic moment of the 
muon (see Section^) only a few perturbative terms are known. Although rapid growth of the perturbative 
coefficients is observed even in relatively low order (see the large number of examples discussed in p8| ), one 
may argue that the factorial growth of the coefficients has not been demonstrated in concrete calculations, 
and that it is unclear if severe cancellations between different classes of diagrams occur in higher order. 

The 30-loop calculation by Broadhurst and Kreimer may indicate in one particular example at least, 
that the factorial divergence is likely to persist, and that possible cancellations due to the renormalization 
or between different sets of diagrams do not contradict the concept of ultimate factorial divergence of the 
perturbative coefficients. From the observation made in [j78| that the S n defined in Eq. (|6^) change little 
at large n, one may tentatively infer the leading factorial divergence of the perturbative coefficients, 

G n ~ 2™- 1 r(n + l/2) n^oo. (66) 

This asymptotic behavior, of course, leads to a vanishing radius of convergence of the perturbative 
expansion (H). 
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Table 3: Perturbative coefficients G n for the anomalous dimension 7 
function of the Yukawa coupling [see Eq. (p^)]. 



n 






_ 

G n 


1 






1 


2 






1 


3 






4 


4 






27 


5 






248 


G 




2 


830 


7 




38 


232 






593 


859 


9 


10 


401 


712 


10 


202 


601 


898 


11 


4 342 


263 


000 


12 


101 551 


822 


350 


13 


2 573 77Q 


506 


192 


14 


70 282 204 


726 


396 


15 


2 057 4Q0 936 


366 


320 


16 


64 291 032 462 


761 


955 


17 


2 136 017 303 903 


513 


184 


18 


75 197 869 250 518 


812 


754 


19 


2 796 475 872 605 709 


079 


512 


20 


1 09 549 71 4 522 464 1 20 


960 


474 


21 


4 509 302 910 783 496 963 


256 


400 


22 


194 584 224 274 515 194 731 


540 


740 


23 


8 784 041 120 771 057 847 338 


352 


720 


24 


414 032 133 398 397 494 698 579 


333 


710 


25 


20 340 342 746 544 244 143 487 152 


873 


888 


2G 


1 039 819 967 521 866 936 447 997 028 


508 


900 


27 


55 230 362 672 853 506 023 203 822 058 


592 


752 


28 


3 043 750 896 574 866 226 650 924 152 479 


935 


036 


29 


173 814 476 864 493 583 374 050 720 641 310 


171 


808 


30 


10 272 611 586 206 353 744 425 870 217 572 111 


879 


288 



As observed by Broadhurst and Kreimer, the perturbation series ([34]) can be resummed with the help 
of an asymptotically improved Borel-Pade technique. From Eq. (22) in |7S[] , it is clear that a (1,1/2)- 
generalized Borel-Pade transformation is used by Broadhurst and Kreimer [for the definition of generalized 
Borel-Pade transformations see Eqs. (||)-(||) in this article]. This is not completely obvious because the 
transformation as used by Broadhurst and Kreimer has been modified additionally such as to normalize 
the first coefficient of the Borel transform (not of the input series) to unity, and the transformation is 
additionally rewritten such as to reflect the vanishing coefficient of zeroth order in a in Eq. (|64|). From 
the leading asymptotics in Eq. (|66|), the singularity of the (1, l/2)-generalized Borel transform closest to 
the origin can be inferred. This singularity was explicitly "put in by hand" by Broadhurst and Kreimer 
(see Eq. (22) in @). 

It is also possible to resum the alternating divergent series (|64|) by a delta transformation, even at 
large coupling. At a large Yukawa coupling of g — 30, we obtain a relative accuracy of 6 significant figures 
in the resummed results with a plain, unmodified delta transformation. 

We add here a remark on the relation of the asymptotically optimized Borel-Pade based predictions 
to those obtained using the delta transformation. We consider the relative accuracy of perturbative pre- 
dictions for the coefficient G30 of the perturbation series defined in Eq. (|64|) using various methods. The 
coefficient G30 is known (see Table |^) , therefore we merely check the accuracy to which this coefficient 
can be reproduced by considering the first 29 perturbative coefficients of the series (p4[). With an asymp- 
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totically optimized (1, l/2)-generalized Borel-Pade technique, the coefficient G 30 can be reproduced with 
a relative accuracy of 5 x 10~ 16 . The delta transformation, without any modifications, leads to a pre- 
diction with a relative error of 3 x 1CT 16 ; this result is more accurate than the prediction provided by 
the asymptotically optimized Borel-Pade transformation. In accordance with the results of Section [|, the 
Borel-Pade transformation can be significantly enhanced by including the pole closest to the origin (see 
Eq. ( |66| ) above and Eq. (22) in j78j). When this pole is included, a prediction is obtained with a relative 
error of 4 x 10~ 17 . We do not consider a-posteriori improvements to either of these predictions, here. 

6 CONCLUSION 

We have considered the resummation of divergent perturbation series and the prediction of unknown 
higher-order perturbative coefficients (perturbative predictions or perturbative extrapolations) . We have 
mentioned and discussed the following resummation prescriptions, 

• the direct application of Pade approximants to a divergent series, 

• the direct application of the nonlinear (delta) sequence transformation, 

• asymptotically improved variants of the Borel-Pade technique. 

The direct application of Pade approximants is less efficient in the resummation of divergent perturbation 
series than both the delta transformation and the combined Borel and Pade techniques. The combined, 
asymptotically improved Borel and Pade techniques, and the delta transformation are complementary. 
On the one hand, it can hardly be overemphasized that the asymptotically improved Borel-Pade tech- 
nique is less general than the delta transformation because it depends on the availability of information 
on the leading large-order asymptotics of the coefficients. By contrast, there is considerable evidence that 
the plain, unmodified delta transformation can sum factorially divergent alternating series which diverge 
as strongly as (3n)\ This is beyond the power of directly applied Pade approximants and also be- 

yond the power of the (1, l)-generalized Borel-Pade transformation ("usual") Borel-Pade transformation 
defined in Eq. (||). 

If additional information is available on the input series, then the asymptotically optimized Borel-Pade 
technique is rather attractive. As discussed in Section |3|, it is possible to enhance the rate of convergence 
simply by utilizing the location of known poles in the Pade approximants to the Borel transform of the 
input series. These improvements are not restricted to the first UV and IR rcnormalon poles, but, as 
shown in Eqs. (|53|)-([55|) and exemplified by the numerical results in Table [j], can take advantage of an 
in principle unlimited number of poles in the Borel plane. Other techniques for possible improvements 
of the Borel-Pade algorithm have been described in |p2| . Note that it is also possible to generalize the 
Borel-Pade technique to those cases where there are poles along the positive real axis in which case the 
Borel integral in Eq. @ is actually undefined J57j. Using the special integration contours in |37|], it is 
even possible to derive nonperturbative imaginary parts from real, not complex, perturbative coefficients. 
Concerning nonperturbative effects in quantum field theory we also refer to the recent investigation |93| . 

As discussed in Section ||, it is possible to accelerate the convergence of the resulting Borel-Pade 
transforms by subsequent application of Wynn's epsilon algorithm (Borel-Pade- Wynn technique). These 
techniques lead to an improved rate of convergence. Note that the use of explicit information of the 
location of the poles in the Borel plane and the subsequent improvement of the convergence of the 
transforms by Wynn's epsilon algorithm should also lead to accelerated convergence in the case of the 
complex integrations discussed in Q . 

Because all the resummation prescriptions discussed above fulfill accuracy-through-order relations, 
they can be used to predict perturbative coefficients (we refer to this procedure as perturbative predic- 
tions or perturbative extrapolations). The straightforward predictions are obtained by reexpansion of the 
rational approximant in powers of the coupling parameter. That is to say, we consider here perturbative 
predictions based on 

• the reexpansion of Pade approximants directly applied to the perturbative (input) series, 

• the reexpansion of nonlinear (delta) sequence transformations directly applied to the perturbative 
(input) series, 
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• and the reexpansion of Pade approximants applied to the asymptotically improved Borel transform 
of the input series. 

As it has been demonstrated in and the predictions based on the delta transformation and on 
the combined Borel and Pade techniques yield better results for the perturbative coefficients of the QED 
effective action than the Pade approximants alone. In jl2| we also presented a number of more realistic 
and practically interesting examples in which the delta transformation leads to better predictions than 
the Pade approximants. 

Note that there is currently no general proof of the assumption that the renormalon poles are the only 
relevant poles in the Borel plane |7^] , but the factorial divergence of the perturbative coefficients is a rather 
commonly accepted assumption |l[, [[7], |8(] We should therefore assume that at least asymptotically, 
the perturbation series in quantum field theory approximate factorially divergent series. This is also 
confirmed by the concrete 30-loop calculation presented in |f78f . For many factorially divergent series, the 
delta transformation produces better numerical results than Pade approximants (see, e.g., Ch. 13 of 0] and 
H |llj). Therefore, the delta transformation can be expected to provide a competitive alternative to Pade 
approximants. We again refer to the large number of recent publications on Pade-based extrapolations 
in quantum field theory (l^-Q . 

We consider here mainly the following asymptotic improvements of perturbative extrapolations 

• a-posteriori corrections based on estimates not only for the next higher-order coefficients, but also 
for the error which is to be expected in the estimation of that coefficient and 



the use of known (renormalon) poles in order to fix the denominatior structure of Pade approximants 
in the context of the Borel-Pade method. 



The general idea of a-posteriori corrections is the following. In higher orders of perturbation theory 
a number of lower-order coefficients are available which may be used in order to construct rational 
approximants. These coefficients can, apart from being useful for the construction the approximants 
itself, also be utilized in order to obtain an estimate for the expected error in the perturbative prediction. 
To this end, the extrapolation procedure is applied to the known lower-order terms in the perturbation 
series. A comparison of the "predictions" for the known lower-order with their exact results gives lower- 
order correction factors which may be extrapolated to higher order. This immediately leads to a correction 
factor for the next higher-order perturbative extrapolation (see Section f|) . 

Using a-posteriori corrections and the pole structure, we obtain improved results for the perturbative 
predictions of the model problem studied in jl2[^5|,|5^] (see Table ||). The a-posteriori corrections also 
lead to an improved estimate for the 10th order anomalous magnetic moment of the muon and bring 
the a-posteriori corrected prediction in close agreement with an analytically obtained estimate [pi| . The 
renormalization-group analysis can lead to a resummation of certain classes of Feynman diagrams, but it 
does not lead to an understanding of topologically new effects which occur in higher orders of perturbation 
theory. This phenomenon has lead to problems in perturbative predictions in the past, especially in those 
cases where these predictions were improved on the basis of a renormalization group analysis. We refer to 
the analysis by Kataev and Starshenko on the 10th order anomalous magnetic moment of the muon ||90f 
and various investigations on the QCD beta function |56||4(|. Notably, as discussed in Section |[ the 
a-posteriori corrected prediction appears to be consistent with the topologically new effects observed in 
10th order of perturbation theory and calculated approximately in ]9l| ]. 

We would like to stress here again that the concept of a-posteriori corrections is rather general and 
can be applied to all prediction algorithms mentioned above. Specifically, we refer to the investigation f[8| 
on significant improvements which can be achieved in the context of Pade-based predictions with this 
technique. The reduction of the magnitude of the a-posteriori correction is an attractive feature of the 
predictions based on the delta transformation. 

In Sections || we show that for the more realistic 30-loop series calculated by Broadhurst and Kreimer 
in J78[ , even the asymptotically optimized Borel-Pade technique cannot quite match the accuracy obtain- 
able by the plain delta transformation. It is only when an additional pole is explicitly put in by hand 
in the Pade transformations that the combined, asymptotically improved Borel-Pade technique becomes 
more accurate than the plain delta transformation. We stress here that the construction of the asymptot- 
ically improved Borel-Pade technique requires in itself a knowledge of the large-order asymptotics of the 
perturbative coefficients. Such information is not available in general cases. Specifically, in those cases 
where only a small number of coefficients are known the leading asymptotics cannot be reliably inferred 
from empirical approaches, either. 
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It has been the purpose of this article to clarify how resummation algorithms and perturbative ex- 
trapolations can be improved if additional information is available on a particular input series. We have 
explained in Sections |^-|5| various algorithms by which the resummation of divergent series and the pre- 
diction of perturbative coefficients can be improved on the basis of additionally available asymptotic 
information on a given series; these improvements can be applied to the Borel-Pade based techniques and 
to the delta transformation techniques [12| . 
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